Comparative genomics of Leishmania donovani progeny from genetic crosses in two sand fly species and impact on the diversity of diagnostic and vaccine candidates.

Sand fly transmitted Leishmania species are responsible for severe, wide ranging, visceral and cutaneous leishmaniases. Genetic exchange can occur among natural Leishmania populations and hybrids can now be produced experimentally, with limitations. Feeding Phlebotomus orientalis or Phlebotomus argentipes on two strains of Leishmania donovani yielded hybrid progeny, selected using double drug resistance and fluorescence markers. Fluorescence activated cell sorting of cultured clones derived from these hybrids indicated diploid progeny. Multilocus sequence typing of the clones showed hybridisation and nuclear heterozygosity, although with inheritance of single haplotypes in a kinetoplastid target. Comparative genomics showed diversity of clonal progeny between single chromosomes, and extraordinary heterozygosity across all 36 chromosomes. Diversity between progeny was seen for the HASPB antigen, which has been noted previously as having implications for design of a therapeutic vaccine. Genomic diversity seen among Leishmania strains and hybrid progeny is of great importance in understanding the epidemiology and control of leishmaniasis. As an outcome of this study we strongly recommend that wider biological archives of different Leishmania species from endemic regions should be established and made available for comparative genomics. However, in parallel, performance of genetic crosses and genomic comparisons should give fundamental insight into the specificity, diversity and limitations of candidate diagnostics, vaccines and drugs, for targeted control of leishmaniasis.

A comparative 'omics' approach for prediction of candidate Strongyloides stercoralis diagnostic coproantigens.

Human infection with the intestinal nematode Strongyloides stercoralis is persistent unless effectively treated, and potentially fatal in immunosuppressed individuals. Epidemiological data are lacking, partially due to inadequate diagnosis. A rapid antigen detection test is a priority for population surveillance, validating cure after treatment, and for screening prior to immunosuppression. We used a targeted analysis of open access 'omics' data sets and used online predictors to identify S. stercoralis proteins that are predicted to be present in infected stool, Strongyloides-specific, and antigenic. Transcriptomic data from gut and non-gut dwelling life cycle stages of S. stercoralis revealed 328 proteins that are differentially expressed. Strongyloides ratti proteomic data for excreted and secreted (E/S) proteins were matched to S. stercoralis, giving 1,057 orthologues. Five parasitism-associated protein families (SCP/TAPS, prolyl oligopeptidase, transthyretin-like, aspartic peptidase, acetylcholinesterase) were compared phylogenetically between S. stercoralis and outgroups, and proteins with least homology to the outgroups were selected. Proteins that overlapped between the transcriptomic and proteomic datasets were analysed by multiple sequence alignment, epitope prediction and 3D structure modelling to reveal S. stercoralis candidate peptide/protein coproantigens. We describe 22 candidates from seven genes, across all five protein families for further investigation as potential S. stercoralis diagnostic coproantigens, identified using open access data and freely-available protein analysis tools. This powerful approach can be applied to many parasitic infections with 'omic' data to accelerate development of specific diagnostic assays for laboratory or point-of-care field application.

Linear and conformational determinants of visceral leishmaniasis diagnostic antigens rK28 and rK39.

BACKGROUND: Recombinant antigens rK39 (based on kinesin sequence) and rK28 (comprising kinesin and HASPB sequences) are a mainstay of serological diagnosis for visceral leishmaniasis (VL). However, their key epitopes and the significance of their structural conformation are not clearly defined, particularly in relation to reported cross-reactivity with sera from patients with malaria, schistosomiasis, and tuberculosis. METHODS: To assess the effect of conformation on antigenicity with Sudanese VL sera, antigens rK39 and rK28 were heat-denatured at 95 °C for 10 min and then assayed by enzyme-linked immunosorbent assay (ELISA). Amino acid sequences of rK39 and rK28 were submitted to NCBI BLASTp to assess homology with Plasmodium, Schistosoma, and Mycobacterium. RESULTS: Heat denaturation significantly diminished the antigenicity of rK39 compared to non-denatured antigen (P = 0.001), but not for rK28 (P = 0.275). In BLASTp searches, HASPB sequences from rK28 had similarities with sequences from Plasmodium, encompassing software-predicted B-cell epitopes. CONCLUSIONS: The antigenicity of rK39 appears to be dependent on structural conformation, whereas that of rK28 depends on linear sequence. HASPB sequence homology with Plasmodium may be responsible for the reported cross-reactivity of rK28 with malaria sera. Further work is warranted to refine the specificity of these antigens.

Microevolution of Trypanosoma cruzi reveals hybridization and clonal mechanisms driving rapid genome diversification.

Protozoa and fungi are known to have extraordinarily diverse mechanisms of genetic exchange. However, the presence and epidemiological relevance of genetic exchange in Trypanosoma cruzi, the agent of Chagas disease, has been controversial and debated for many years. Field studies have identified both predominantly clonal and sexually recombining natural populations. Two of six natural T. cruzi lineages (TcV and TcVI) show hybrid mosaicism, using analysis of single-gene locus markers. The formation of hybrid strains in vitro has been achieved and this provides a framework to study the mechanisms and adaptive significance of genetic exchange. Using whole genome sequencing of a set of experimental hybrids strains, we have confirmed that hybrid formation initially results in tetraploid parasites. The hybrid progeny showed novel mutations that were not attributable to either (diploid) parent showing an increase in amino acid changes. In long-term culture, up to 800 generations, there was a variable but gradual erosion of progeny genomes towards triploidy, yet retention of elevated copy number was observed at several core housekeeping loci. Our findings indicate hybrid formation by fusion of diploid T. cruzi, followed by sporadic genome erosion, but with substantial potential for adaptive evolution, as has been described as a genetic feature of other organisms, such as some fungi.

Congenital Chagas disease in Santa Cruz Department, Bolivia, is dominated by Trypanosoma cruzi lineage V.

BACKGROUND: This study identified Trypanosoma cruzi discrete typing units (DTUs) in maternal and infant specimens collected from two hospitals in Bolivia, using conventional genotyping and DTU-specific serotyping. METHODS: Specimens from 142 mothers were used, including 24 seronegative and 118 seropositive individuals; 29 women transmitted T. cruzi to their infants. Maternal and infant parasite loads were determined by quantitative real-time PCR. Maternal sera were tested with an in-house parasite lysate ELISA and serotyped by a lineage-specific peptide ELISA, targeting the trypomastigote small surface antigen (TSSA). Trypanosoma cruzi genotypes in infected infants were determined by a triple PCR-RFLP assay. RESULTS: All infant specimens were genotyped as TcV. Maternal parasite loads and absorbance values by the lysate ELISA were significantly higher for transmitters compared with non-transmitters. Among seropositive mothers, 65.3% had positive results by the TSSA II/V/VI peptide ELISA. No significant difference in reactivity to TSSA II/V/VI was observed for transmitters compared with non-transmitters (79.3% vs 60.7%, respectively). CONCLUSIONS: Our findings reinforce the difficulty in obtaining sufficient sample numbers and parasite DNA to investigate the interaction between parasite genetics and the risk of congenital transmission and argue for the inclusion of DTU-specific serotyping in prospective studies.

Assessing antibody decline after chemotherapy of early chronic Chagas disease patients.

BACKGROUND: Chagas disease remains a significant public health problem in Latin America. There are only two chemotherapy drugs, nifurtimox and benznidazole, and both may have severe side effects. After complete chemotherapy of acute cases, seropositive diagnosis may revert to negative. However, there are no definitive parasitological or serological biomarkers of cure. METHODS: Following a pilot study with seven Bolivian migrants to Spain, we tested 71 serum samples from chronic patients (mean age 12.6 years) inhabiting the Argentine Chaco region. Benznidazole chemotherapy (5-8 mg/kg day, twice daily for 60 days) was administered during 2011-2016. Subsequently, pre-and post-chemotherapy serum samples were analysed in pairs by IgG1 and IgG ELISA using two different antigens and Chagas Sero K-SeT rapid diagnostic tests (RDT). Molecular diagnosis by kDNA-PCR was applied to post-treatment samples. RESULTS: Pilot data demonstrated IgG1 antibody decline in three of seven patients from Bolivia 1 year post-treatment. All Argentine patients in 2017 (averaging 5 years post-treatment), except one, were positive by conventional serology. All were kDNA-PCR-negative. Most (91.5%) pre-treatment samples were positive by the Chagas Sero K-SeT RDT, confirming the predominance of TcII/V/VI. IgG1 and IgG of Argentine patients showed significant decline in antibody titres post-chemotherapy, with either lysate (IgG, P = 0.0001, IgG1, P = 0.0001) or TcII/V/VI peptide antigen (IgG, P = 0.0001, IgG1, P = 0.0001). IgG1 decline was more discriminative than IgG. Antibody decline after treatment was also detected by the RDT. Incomplete treatment was associated with high IgG1 post-treatment titres against lysate (P = 0.013), as were IgG post-treatment titres to TcII/V/VI peptide (P = 0.0001). High pre-treatment IgG1 with lysate was associated with Qom ethnicity (P = 0.045). No associations were found between gender, age, body mass index and pre- or post-treatment antibody titres. CONCLUSIONS: We show that following chemotherapy of early chronic Chagas disease, significant decline in IgG1 antibody suggests cure, whereas sustained or increased IgG1 is a potential indicator of treatment failure. Due to restricted sensitivity, IgG1 should not be used as a diagnostic marker but has promise, with further development, as a biomarker of cure. We show that following chemotherapy of early chronic Chagas disease, a significant decline in IgG1 antibody suggests cure, whereas sustained or increased IgG1 is a potential indicator of treatment failure. Due to restricted sensitivity, IgG1 should not be used as a diagnostic marker but has promise, with further development, as a biomarker of cure.

Emergent canine visceral leishmaniasis in Argentina: Comparative diagnostics and relevance to proliferation of human disease.

BACKGROUND: Visceral leishmaniasis (VL) is a zoonotic protozoal vector-borne disease that is a major public health challenge. In Argentina, canine (CVL) and human visceral leishmaniasis (HVL) have recently emerged. There is a lack of standardised diagnostic tests for CVL, which hinders control of CVL and HVL. METHODOLOGY/PRINCIPAL FINDINGS: Sampling was carried out in Puerto Iguazú, Argentina, comprising 190 asymptomatic, oligosymptomatic and polysymptomatic dogs. The following diagnostics were applied: microscopy of lymph node aspirate (LNA); three immunochromatographic rapid diagnostic tests (RDTs), prototype rK28-ICT, rK39-ICT (both Coris BioConcept), commercial rK39 (InBios); ELISA for IgG, IgG1 and IgG2, against rK28, rK39 or crude lysate antigen. DNA detection and analysis, with 30 dogs, was of the ITS1 region using skin samples, and loop-mediated isothermal amplification (LAMP; Eiken Loopamp) of buffy coat, skin scrape or LNA. 15.4% of dogs were positive by LNA microscopy. The rK28 RDT had higher seropositivity rate (61%) than either a prototype rK39 RDT (31.4%) or commercial rK39 RDT (18.8%), without cross-reactivity with six other pathogens. IgG anti-rK39 ELISA antibody titres, but not IgG2, were positively correlated with number of clinical signs. LAMP with LNA had a higher positivity rate than PCR; buffy coat sampling was more sensitive than skin scrape. ITS1 confirmed Leishmania (Leishmania) infantum as the agent of CVL. Leishmania (Viannia) spp. was detected in skin samples from two dogs, compatible with Leishmania (Viannia) braziliensis. CONCLUSIONS/SIGNIFICANCE: Seroprevalence confirmed rapid increase in CVL in Puerto Iguazú. The rK28 RDT test potentially has great value for improved point-of-care diagnosis. Given cost reduction and accessibility, commercial LAMP may be applicable to buffy coat. RDT biomarkers of CVL clinical status are required to combat spread of CVL and HVL. The presence of Viannia, perhaps as an agent of human mucocutaneous leishmaniasis (MCL), highlights the need for vigilance and surveillance.

Under pressure: phenotypic divergence and convergence associated with microhabitat adaptations in Triatominae.

BACKGROUND: Triatomine bugs, the vectors of Chagas disease, associate with vertebrate hosts in highly diverse ecotopes. It has been proposed that occupation of new microhabitats may trigger selection for distinct phenotypic variants in these blood-sucking bugs. Although understanding phenotypic variation is key to the study of adaptive evolution and central to phenotype-based taxonomy, the drivers of phenotypic change and diversity in triatomines remain poorly understood. METHODS/RESULTS: We combined a detailed phenotypic appraisal (including morphology and morphometrics) with mitochondrial cytb and nuclear ITS2 DNA sequence analyses to study Rhodnius ecuadoriensis populations from across the species' range. We found three major, naked-eye phenotypic variants. Southern-Andean bugs primarily from vertebrate-nest microhabitats (Ecuador/Peru) are typical, light-colored, small bugs with short heads/wings. Northern-Andean bugs from wet-forest palms (Ecuador) are dark, large bugs with long heads/wings. Finally, northern-lowland bugs primarily from dry-forest palms (Ecuador) are light-colored and medium-sized. Wing and (size-free) head shapes are similar across Ecuadorian populations, regardless of habitat or phenotype, but distinct in Peruvian bugs. Bayesian phylogenetic and multispecies-coalescent DNA sequence analyses strongly suggest that Ecuadorian and Peruvian populations are two independently evolving lineages, with little within-lineage phylogeographic structuring or differentiation. CONCLUSIONS: We report sharp naked-eye phenotypic divergence of genetically similar Ecuadorian R. ecuadoriensis (nest-dwelling southern-Andean vs palm-dwelling northern bugs; and palm-dwelling Andean vs lowland), and sharp naked-eye phenotypic similarity of typical, yet genetically distinct, southern-Andean bugs primarily from vertebrate-nest (but not palm) microhabitats. This remarkable phenotypic diversity within a single nominal species likely stems from microhabitat adaptations possibly involving predator-driven selection (yielding substrate-matching camouflage coloration) and a shift from palm-crown to vertebrate-nest microhabitats (yielding smaller bodies and shorter and stouter heads). These findings shed new light on the origins of phenotypic diversity in triatomines, warn against excess reliance on phenotype-based triatomine-bug taxonomy, and confirm the Triatominae as an informative model system for the study of phenotypic change under ecological pressure .

Repeat-Driven Generation of Antigenic Diversity in a Major Human Pathogen, Trypanosoma cruzi.

Trypanosoma cruzi, a zoonotic kinetoplastid protozoan parasite, is the causative agent of American trypanosomiasis (Chagas disease). Having a very plastic, repetitive and complex genome, the parasite displays a highly diverse repertoire of surface molecules, with pivotal roles in cell invasion, immune evasion and pathogenesis. Before 2016, the complexity of the genomic regions containing these genes impaired the assembly of a genome at chromosomal level, making it impossible to study the structure and function of the several thousand repetitive genes encoding the surface molecules of the parasite. We here describe the genome assembly of the Sylvio X10/1 genome sequence, which since 2016 has been used as a reference genome sequence for T. cruzi clade I (TcI), produced using high coverage PacBio single-molecule sequencing. It was used to analyze deep Illumina sequence data from 34 T. cruzi TcI isolates and clones from different geographic locations, sample sources and clinical outcomes. Resolution of the surface molecule gene distribution showed the unusual duality in the organization of the parasite genome, a synteny of the core genomic region with related protozoa flanked by unique and highly plastic multigene family clusters encoding surface antigens. The presence of abundant interspersed retrotransposons in these multigene family clusters suggests that these elements are involved in a recombination mechanism for the generation of antigenic variation and evasion of the host immune response on these TcI strains. The comparative genomic analysis of the cohort of TcI strains revealed multiple cases of such recombination events involving surface molecule genes and has provided new insights into T. cruzi population structure.

Glycosylation of Trypanosoma cruzi TcI antigen reveals recognition by chagasic sera.

Chagas disease is considered the most important parasitic disease in Latin America. The protozoan agent, Trypanosoma cruzi, comprises six genetic lineages, TcI-TcVI. Genotyping to link lineage(s) to severity of cardiomyopathy and gastrointestinal pathology is impeded by the sequestration and replication of T. cruzi in host tissues. We describe serology specific for TcI, the predominant lineage north of the Amazon, based on expression of recombinant trypomastigote small surface antigen (gTSSA-I) in the eukaryote Leishmania tarentolae, to allow realistic glycosylation and structure of the antigen. Sera from TcI-endemic regions recognised gTSSA-I (74/146; 50.7%), with no cross reaction with common components of gTSSA-II/V/VI recombinant antigen. Antigenicity was abolished by chemical (periodate) oxidation of gTSSA-I glycosylation but retained after heat-denaturation of conformation. Conversely, non-specific recognition of gTSSA-I by non-endemic malaria sera was abolished by heat-denaturation. TcI-specific serology facilitates investigation between lineage and diverse clinical presentations. Glycosylation cannot be ignored in the search for immunogenic antigens.

Central Asian Rodents as Model Animals for Leishmania major and Leishmania donovani Research.

The clinical manifestation of leishmaniases depends on parasite species, host genetic background, and immune response. Manifestations of human leishmaniases are highly variable, ranging from self-healing skin lesions to fatal visceral disease. The scope of standard model hosts is insufficient to mimic well the wide disease spectrum, which compels the introduction of new model animals for leishmaniasis research. In this article, we study the susceptibility of three Asian rodent species (Cricetulus griseus, Lagurus lagurus, and Phodopus sungorus) to Leishmania major and L. donovani. The external manifestation of the disease, distribution, as well as load of parasites and infectiousness to natural sand fly vectors, were compared with standard models, BALB/c mice and Mesocricetus auratus. No significant differences were found in disease outcomes in animals inoculated with sand fly- or culture-derived parasites. All Asian rodent species were highly susceptible to L. major. Phodopus sungorus showed the non-healing phenotype with the progressive growth of ulcerative lesions and massive parasite loads. Lagurus lagurus and C. griseus represented the healing phenotype, the latter with high infectiousness to vectors, mimicking best the character of natural reservoir hosts. Both, L. lagurus and C. griseus were also highly susceptible to L. donovani, having wider parasite distribution and higher parasite loads and infectiousness than standard model animals.

Isolation and characterisation of Leishmania donovani protein antigens from urine of visceral leishmaniasis patients.

Diagnosis of visceral leishmaniasis (VL) relies on invasive and risky aspirate procedures, and confirmation of cure after treatment is unreliable. Detection of Leishmania donovani antigens in urine has the potential to provide both a non-invasive diagnostic and a test of cure. We searched for L. donovani antigens in urine of VL patients from India and Sudan to contribute to the development of urine antigen capture immunoassays. VL urine samples were incubated with immobilised anti-L. donovani polyclonal antibodies and captured material was eluted. Sudanese eluted material and concentrated VL urine were analysed by western blot. Immunocaptured and immunoreactive material from Indian and Sudanese urine was submitted to mass spectrometry for protein identification. We identified six L. donovani proteins from VL urine. Named proteins were 40S ribosomal protein S9, kinases, and others were hypothetical. Thirty-three epitope regions were predicted with high specificity in the 6 proteins. Of these, 20 were highly specific to Leishmania spp. and are highly suitable for raising antibodies for the subsequent development of an antigen capture assay. We present all the identified proteins and analysed epitope regions in full so that they may contribute to the development of non-invasive immunoassays for this deadly disease.

Novel structural CYP51 mutation in Trypanosoma cruzi associated with multidrug resistance to CYP51 inhibitors and reduced infectivity.

Ergosterol biosynthesis inhibitors, such as posaconazole and ravuconazole, have been proposed as drug candidates for Chagas disease, a neglected infectious tropical disease caused by the protozoan parasite Trypanosoma cruzi. To understand better the mechanism of action and resistance to these inhibitors, a clone of the T. cruzi Y strain was cultured under intermittent and increasing concentrations of ravuconazole until phenotypic stability was achieved. The ravuconazole-selected clone exhibited loss in fitness in vitro when compared to the wild-type parental clone, as observed in reduced invasion capacity and slowed population growth in both mammalian and insect stages of the parasite. In drug activity assays, the resistant clone was above 300-fold more tolerant to ravuconazole than the sensitive parental clone, when the half-maximum effective concentration (EC50) was considered. The resistant clones also showed reduced virulence in vivo, when compared to parental sensitive clones. Cross-resistance to posaconazole and other CYP51 inhibitors, but not to other antichagasic drugs that act independently of CYP51, such as benznidazole and nifurtimox, was also observed. A novel amino acid residue change, T297M, was found in the TcCYP51 gene in the resistant but not in the sensitive clones. The structural effects of the T297M, and of the previously described P355S residue changes, were modelled to understand their impact on interaction with CYP51 inhibitors.

Application of a recombinase polymerase amplification (RPA) assay and pilot field testing for Giardia duodenalis at Lake Albert, Uganda.

BACKGROUND: Giardia duodenalis is a gastrointestinal protozoan causing 184 million cases of giardiasis worldwide annually. Detection is by microscopy or coproantigen assays, although sensitivity is often compromised by intermittent shedding of cysts or trophozoites, or operator expertise. Therefore, for enhanced surveillance field-applicable, point-of-care (POC), molecular assays are needed. Our aims were to: (i) optimise the recombinase polymerase amplification (RPA) assay for the isothermal amplification of the G. duodenalis ß-giardin gene from trophozoites and cysts, using published primer and probes; and (ii) perform a pilot field validation of RPA at a field station in a resource-poor setting, on DNA extracted from stool samples from schoolchildren in villages around Lake Albert, Uganda. Results were compared to an established laboratory small subunit ribosomal RNA (SSU rDNA) qPCR assay with additional testing using a qPCR targeting the triose phosphate isomerase (tpi) DNA regions that can distinguish G. duodenalis of two different assemblages (A and B), which are human-specific. RESULTS: Initial optimisation resulted in the successful amplification of predicted RPA products from G. duodenalis-purified gDNA, producing a double-labelled amplicon detected using lateral flow strips. In the field setting, of 129 stool samples, 49 (37.9%) were positive using the Giardia/Cryptosporidium QuikChek coproantigen test; however, the RPA assay when conducted in the field was positive for a single stool sample. Subsequent molecular screening in the laboratory on a subset (n = 73) of the samples demonstrated better results with 21 (28.8%) RPA positive. The SSU rDNA qPCR assay resulted in 30/129 (23.3%) positive samples; 18 out of 73 (24.7%) were assemblage typed (9 assemblage A; 5 assemblage B; and 4 mixed A+B). Compared with the SSU rDNA qPCR, QuikChek was more sensitive than RPA (85.7 vs 61.9%), but with similar specificities (80.8 vs 84.6%). In comparison to QuikChek, RPA had 46.4% sensitivity and 82.2% specificity. CONCLUSIONS: To the best of our knowledge, this is the first in-field and comparative laboratory validation of RPA for giardiasis in low resource settings. Further refinement and technology transfer, specifically in relation to stool sample preparation, will be needed to implement this assay in the field, which could assist better detection of asymptomatic Giardia infections.

Genomic analysis of natural intra-specific hybrids among Ethiopian isolates of Leishmania donovani.

Parasites of the genus Leishmania (Kinetoplastida: Trypanosomatidae) cause widespread and devastating human diseases. Visceral leishmaniasis due to Leishmania donovani is endemic in Ethiopia where it has also been responsible for major epidemics. The presence of hybrid genotypes has been widely reported in surveys of natural populations, genetic variation reported in a number of Leishmania species, and the extant capacity for genetic exchange demonstrated in laboratory experiments. However, patterns of recombination and the evolutionary history of admixture that produced these hybrid populations remain unclear. Here, we use whole-genome sequence data to investigate Ethiopian L. donovani isolates previously characterized as hybrids by microsatellite and multi-locus sequencing. To date there is only one previous study on a natural population of Leishmania hybrids based on whole-genome sequences. We propose that these hybrids originate from recombination between two different lineages of Ethiopian L. donovani occurring in the same region. Patterns of inheritance are more complex than previously reported with multiple, apparently independent, origins from similar parents that include backcrossing with parental types. Analysis indicates that hybrids are representative of at least three different histories. Furthermore, isolates were highly polysomic at the level of chromosomes with differences between parasites recovered from a recrudescent infection from a previously treated individual. The results demonstrate that recombination is a significant feature of natural populations and contributes to the growing body of data that shows how recombination, and gene flow, shape natural populations of Leishmania.

Correction: A lineage-specific rapid diagnostic test (Chagas Sero K-SeT) identifies Brazilian Trypanosoma cruzi II/V/VI reservoir hosts among diverse mammalian orders.

[This corrects the article DOI: 10.1371/journal.pone.0227828.].

Global genome diversity of the Leishmania donovani complex.

Protozoan parasites of the Leishmania donovani complex - L. donovani and L. infantum - cause the fatal disease visceral leishmaniasis. We present the first comprehensive genome-wide global study, with 151 cultured field isolates representing most of the geographical distribution. L. donovani isolates separated into five groups that largely coincide with geographical origin but vary greatly in diversity. In contrast, the majority of L. infantum samples fell into one globally-distributed group with little diversity. This picture is complicated by several hybrid lineages. Identified genetic groups vary in heterozygosity and levels of linkage, suggesting different recombination histories. We characterise chromosome-specific patterns of aneuploidy and identified extensive structural variation, including known and suspected drug resistance loci. This study reveals greater genetic diversity than suggested by geographically-focused studies, provides a resource of genomic variation for future work and sets the scene for a new understanding of the evolution and genetics of the Leishmania donovani complex.